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Abstract 

In this chapter we discuss the X-ray radiation from relativistic accretion disks 
around supermassive black holes, supposed to exist in the centers of Active 
Galactic Nuclei (AGN). Our focus is on the X-ray radiation, especially in 
the Fe Ka line which originates in the innermost parts of an accretion disk. 
Moreover, here we discuss some effects which can disturb the Fe Ka profile 
and cause its rapid and irregular variability, observed in the X-ray spectra 
of some AGN. We will pay attention to three such effects: perturbations 
in the disk emissivity, absorbtion by warm absorbers and gravitational mi- 
crolensing. The X-ray emission from accretion disks around non-rotating 
(Schwarzschild metric), as well as rotating (Kerr metric) supermassive black 
holes, is discussed. The X-ray radiation of AGN is probably produced in a 
compact region near their central supermassive black holes, and can provide 
us some essential information about the plasma conditions and the space- 
time geometry in these regions. The goal of this chapter is mainly to present 
a short overview of some important and recent investigations in this field. 



1 Introduction 



Active Galactic Nuclei are powerful sources of radiation in a wide spectral 
range: from 7 rays to radio waves. They derive their extraordinary lu- 
minosities (sometimes more than 10.000 times higher than luminosities of 
"ordinary" galaxies) from energy release by matter a ccreting towards , and 
falling into, a central supermassive black hole ( see e.g. iPetersoru, |2004| ). 

According to the unification model of AGN (lAntonuccil . Il993l ). their cen- 
tral engine consists of a supermass ive black hole (SMBH) with mass rangiri g 



from 10^ to 10^ solar masses Mq flKaspi et all . I2OOOI : IPeterson et all . I2OOJ), 



which is surrounded by an accretion disk that radiates in the X-ray band 
(IFabian et al.l . 119891 ). The vast majority of the X-ray sources in the Universe 
are AGN. The integral emission of AGN reflects the history of accretion onto 
SMBH over cosmic time. Emission lines are usually seen in the X-ray spec- 
tra of AGN. The broad emission Fe Ka spectral line (6.4-6.9 keV, depending 
on ionization state) with asymmetric profile (narrow bright blue peak and 
a w ide faint red wing) h as been observed in a number of type 1 AGN (see 



e.g. iNandra et al. . 20071 ). Early results from the ASCA (Advanced Satellite 
for Cosmology and Astrophysics) era suggested that broad relativistic lines 
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might be common in type 1 AGN. However, surprisingly, they have been 
detected in significant nu mber and their features d escribed in only a small 
fraction of those sources (IReynolds &: Nowakl. 120031) . The first and the bes t 
studied one is MCG-6-30-15 flXanaka et all Il995l : iFabian fc Vaughanl . l2003h . 

In some cases the line width corresponds to one third of speed of light, 
indicating that its emitters rotate with relativistic velocities. Therefore, the 
line is probably produced in a very compact region near the central black 
hole of AGN and can provide us some essential information about the plasma 
conditions and the space-time geometry in vicinity of the black hole. Conse- 
quently, if the line is emitted close enou gh to the SMBH, it shows a broad rel- 
ativistic profile affected by SMBH spin (IReynolds fc Fabianl. [ioosh and grav- 



itational redshift (as well as other general relativistic effects) ( iFabian et al. 



19891 : iLaod . llQQlD 



Several studies have been performed over samples of loca l AGN (see e.g . 



Nandra et all. 119971: lYaaoob et all . l2005l : iNandra et al.l . 120071 : iBianchi et al 



200q : lMarkowitz et al.l . l2008l . etc), as well as from di stant quasars fICorral et al. 



20081 ) in order to characterize the Fe Ka emission. iNandra et al.l (120071 ) per- 
formed a spectral analysis of a sample of 26 type 1 to 1.9 Seyferts galaxies 
{z < 0.05) observed by XMM-Newton. They found that a relativistic line is 
significantly detected in a half of their sample (54±10 percent) with a mean 
equivalent width (EW) of ~ 80 eV, but around 30% of selected AGN showed 
a relativistic broad line that can be explained by the emission of an accretion 
disk. 

Accretion disks could have different forms, dimensions, and emission, de- 
pending on the type of central black hole, whether it is rotating (Kerr metric) 
or non-rotating (Schwarzschild metric). They represent an efficient mecha- 
nism for extracting gravitational potential energy and converting it into ra- 
diation, giving us the most probable explanation for the main characteristics 
of AGN (high luminosity, compactness, jet formation, rapid time variations 
in radiation and the profile of the Fe Ka spectral line). 

Here we discuss the X-ray radiation from relativistic accretion disks around 
supermassive black holes, supposed to exist in the centers of AGN. Especially, 
we discuss the Fe Ka line profile which originates from the accretion disk. 
Moreover, we also present some effects which can disturb the Fe Ka profile, 
such as: perturbations in the disk emissivity, absorbtion by warm absorbers 
and gravitational microlensing. 

The aim of this chapter is to present a short overview of results of some 
recent investigations in this field, and it is divided into the following six 
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sections: Active Galactic Nuclei as Hosts of Supermassive Black Holes - 
where the main features, classification and unified model of active galax- 
ies are briefly presented, Space-Time Geometry in Vicinity of Supermassive 
Black Holes - where the basic definitions of Schwarzschild and Kerr met- 
rics are given. Accretion Disk Around a Supermassive Black Hole - where 
we explain the standard model of an accretion disk, including its emission, 
accretion rate, luminosity, structure and spectral distribution, Supermassive 
Black Holes and X-ray Emission - where the focus is on the modeling of 
the observed X-ray radiation from relativistic accretion disk around a su- 
permassive black hole of AGN in both the Fe Ka spectral line and X-ray 
continuum. Variability of X-ray Emission Around Supermassive Black Hole 
- where we present the main causes of rapid and irregular variability of the 
X-ray emission which can be due to disk instability, reflecting in perturba- 
tions of its emissivity, or it could be caused by some external effects, such 
as gravitational microlensing and absorption by X-ray absorbers, and finally 
Gonclusion - where the most important results from previous sections are 
pointed out and their brief summary is given. 

Finally, we should note that a huge number of papers devoted to investi- 
gation of the X-ray emission of AGN was published in the last two decades, 
so it was not possible to mention all of them here, but more details and refer - 
ences c an be f ound in review papers s uch a s e.g. iBrandt fc Hasingen (120051 ): 



Millerl (120071 1: Harris fc Krawczvnskil (12006[ ). etc. 



2 Active Galactic Nuclei as Hosts of Super- 
massive Black Holes 

Active galaxies differ from the so called "normal" galaxies in the amount of 
energy emitted from their nuclei. The term Active Galactic Nuclei (AGN) 
refers to the energetic phenomena in the central region of a galaxy which 
cannot be solely or directly produced by stars. The nebular-like emission 
spectra of NGC 1068, NGC 4051 and NGC 4151, with Ha, [ O II] A3727, [N e 



III] A3869 and [O III] AA 4363,4959,5007 were observed bv iHubbld (119261) . 



but the first classification of these objects was made by ISeyfertI (119431 ). He 
recognized these objects as a class of galaxies with strong and high-excitation 
optical emission lines localized in the nuclei^. 



^For more detailed historical review see 



Osterbrockl (Il989f l 
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Active galaxies have been the subject of an intensive astrophysical in- 
vestigations for the last 3-4 decades. As the most luminous objects with 
one of the most powerful energy release rates and with the most compact 
dimensions, AGN are of the great interest in modern astrophysics. AGN (i.e. 
quasars) are the brightest objects in the Universe and thus, these objects 
are important for the studies of the early Universe and cosmology in general. 
For instance, studies of their luminosity functions or quasar host galaxies are 
crucial for understanding the formation and evolution of galaxies in general. 
The class of active galaxies contains many different objects, such as quasars, 
gigantic elliptic radio galaxies, luminous spiral Seyfert galaxies, blazars, etc. 



2.1 Main characteristics of AGN 



AGN are the most luminous objects in the Universe and they have luminosi- 
ties in the range from ~ 10^^ to 10^^ erg s~^, that can be 10^ times greater 
than in the case of a typical galaxy. Apart from their great luminosities, 
AGN are emitting in the broad band of electromagnetic spectrum. The ob- 
served emission of AGN is in the continuum and in l ines, from the 7 and X 
domains to the far infrared and radio bands 



PetersonI, 12Q0J) 



[see e.g. 

One of the characteristics of AGN is violent and fast variability observed 
in different parts of electromagnetic spectra. Fast changes in the AGN bright- 
ness are observed s ometimes just for a coup le of days, e.g. in the case of 



galaxy NGC 4151 (iShapovalova et al.l . |2008| ). This leads to the conclusion 



that the emission regions in the AGN have small dimensions (from only cou- 
ple of light days up to couple of light months). Also, the brightness of AGN 
can be so high that the brightness of the host galaxy can be neglected in the 
total brightness (for instance in the case of quasars). Consequently, they are 
compact regions which could not be directly observed in the most bands of 
electromagnetic spectrum. Therefore, AGN and their physical/kinematical 
properties are indirectly studied through the analysis of their spectra. 

In contrast to the spectral energy distribution (SED) of an ordinary 
galaxy (that represents a sum of stellar spectra, thus the most of its lu- 
minosity comes within no more than one decade of frequency), AGN emit 
across the broad range of frequencies (of order 10^ Hz). The AGN contin- 
uum in a number of high-luminosity s ources peaks in the ultraviolet, i.e. in 
the so called Big Blue Bump (BBB) fIShildsl . Il978l : Izhou et all . Il997h . The 



BBB is attributed to the thermal emission from an optically thick region of 
the accretion disk fIShildsl . Il978l : luirich et al.1 . Il980h . or from optically thin 
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regio ns, i.e. free-free emission (lAntonucci fc Barvainid . Il988l : iFerland et al. 
The observed SED in AGN is remarkably diff erent from the therma l 



(black-body) spectrum of a star or a regular galaxy (lOke fc Sargentl . Il968l ). 
This featureless continuum observed in AGN characterizes the main source 
of energy input and suggests that the mechanisms that produce it are com- 
mon to all types of AGN. The SED can be approximated with the power 
l aw f unction F where a is the spectral index (iKrolikl . Il999l : iPetersoru . 



20041 ). 



The observed continuum of a typical Seyfert 1 galaxy is very strong com- 
pared with Seyfert 2 galaxies; as a result Seyfert 1 galaxies appear to be more 
luminous than Seyfert 2-s. Some AGN have the SED almost flat from the 
infrared to X-ray part of the spectrum, so that the spectral index is a ~ 1, 
although it is usually steeper. In the radio band, the brightness of active 
galaxies is usually an order of magnitude higher than of normal galaxies. 
But still, even if it is so intensive, the radio luminosity is never higher more 
than 1% of the bolometric luminosity. 

Another property of AGN is strong X-ray flux, with the X-ray continuum 
produced by lower energy photo ns which are Compton scatt ered to higher 
energies by relativistic electrons (ISunyaev &: Titarchukl . Il980l ). The fraction 
of the power emitted in the X- ray ena i ssion i s three o r four times larger in 
AGN than in normal galaxies (iKrolikl . Il999l : iPetersonl . 120041 ). therefore the 
X-emission can indicate the presence of AGN. 

Today, it is widely accepted that the mechanism which powers AGN is 
accretion onto supermassive black hole (SMBH), i.e. the greatest part of 
the continuum emission comes from an accretion disk and its corona. The 
mass of the supermas s ive black hole is estimated to be from 10^ to 10^ Mq 



flPeterson et all . l2004l : IPeterson fc Wandej . l2000h . Thus, in the center of an 



AGN there is a supermassive black hole with an accretion disk surrounded 
by gas and dust in the form of a torus. 

There are also other emitting regions present in the vicinity of the super- 
massive black hole, such as ionized gas clouds that produce intensive emission 
lines or jets of matter that are mostly visible in the radio band, but some- 
times they can also be seen in the optical band. The kinematics of these 
regions is very complex as well as the physical processes that produce such 
specific spectra of these objects. In the central part, there are mainly three 
emitting line regions: Fe Ka (in the X-ray range j^. Broad and Narrow Line 



^See §5.2 for more details 
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Emitting Regions (BLR and NLR, respectively; in optical/UV spectra). 

In the optical and UV band, the total flux of emission (and occasionally 
absorbtion) lines contributes from several percent to tens of percent of the 
continuum flux. The existence of the broad and narrow emission lines comes 
from the fact that an AGN contains two separated emission line regions 
with different kinematics, density, ionization, optical thickness and radiation 
transfer. 

The basic characteristics of the BLR and NLR can be described as follows: 



BLR is a compact region, with di mensions from only couple of light 



days up to couple of light weeks (iKaspi et al.l . l2000l : iPeterson et al 



20041 ). that is located in the vicinity of the black hole (<1 pc). The 
structure of this region is very complex and most likely it c onsist s 
from at l e ast tw o separated subreg i ons (see e.g Popovic et al.l (2004); 
Ihc et all fl2006h : iBon et al.l fl2006h : IPopovic et al.l fl2008h . etc). The 
ionized gas in this region is of relatively high density {n > lO^cm^^) 
and temperature T ~ 10^ K. Emitters are moving with high velocities, 
up to 10000 km/s. In this region the broad emission lines (BELs) from 
allowed transitions are formed. 



NLR extends even up to 1 kpc far from the central source f Petersonl . 
20041 ). The density in this region is significantly smaller (n ~ lO^cm"^) 
than in the BLR, while the temperature is on the same order of magni- 
tude as in the BLR. From this region the narrow emission lines (NELs) 
arise, indicating a random motion of the emitting gas smaller than 
1000 kms~^. NELs are often coming from the forbidden transitions. 
The physical properties of this region are much closer to the properties 
of emission nebulae than to the BLR. 



According to the different characteristics mentioned above, there are sev- 
eral classes of AGN. 



2.2 Classification 

The first group of stellar systems with active nuclei consists of Seyfert galax- 
ies. These are spiral galaxies with very bright nuclei which spectra have 
strong emission lines of neutral and multi-ionized emitters. The presence of 
highly ionized emission lines indicates the existence of a non-stellar ionization 
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continuum. The presence or absence of BELs (i.e. BLR) has been histor- 
ically used to separate Seyfert galaxies into two classes. Seyfert 1 galaxies 
have broad permitted (H I, He I and He II) and narrow permitted and forbid- 
den lines (such as [O HI] A500 7), while Seyfert 2 galaxi e s hav e only narrow 
permitted and forbidden lines (IKhachikian fc Weedmanl . Il974l ) . This simple 
Seyfert classification scheme can be fu rther sub-divided according to specific 
spectral properties (jOsterbroda . Il989l ). as e.g. Seyfert galaxies with inter- 
mediate Balmer profiles are classified as Seyfert 1.5-s with apparent narrow 
Balmer line components superimposed on broad wings. 

Many low luminosity galaxies have a nucleus that resembles the Seyfert 
2 nucleus, but with the spectrum that shows the forbidden lines from lower 
ionization states. Therefore, these galaxies are called Low-Ionization Nuclear 
Emission-line Regions, i.e LINERs. 

Active galaxies which emit intensively in the radio band are so called radio 
galaxies (RG). Radio galaxies are dominantly elliptic galaxies that are divided 
into two morphologically different types: PR I and PR II. The classification 
is done according to the relative position of the surface brightness maximum 
(called the hot spot) in the radio lobes. Por example, in the PR I radio galaxy 
the distance between the hot spot in the radio lobes is less than half of the 
maxim um diameter of the radio sou rce, while in the PR II type this distance is 
larger (IBinney fc Merrifieldl . Il998l ). Like Seyfert galaxies, radio galaxies can 
be divided into subgroups according to the widths of their emission hues, 
i.e. Narrow Line Radio Galaxies (NLRG) that emit narrow emission lines 
characteristic for Seyfert 2 galaxies, and Broad Line Radio Galaxies (BLRG) 
that emit broad lines as Seyfert 1 galaxies. 

Being the most powerful sources in the Universe, quasars b elong to the 
group of active galaxies with the absolute magnitude Mb < —23 (lOsterbrockl . 
mm. These quasi stellar (quasars) objects emit intensively in the radio do- 



main and their spectrum shows broad emission lines significantly redshifted. 
In this group of AGN there are also Quasi-Stellar Objects (QSOs) which 
have the same observed properties as quasars, except they are weak radio 
sources. Besides, there is a group of QSOs with broad absorption lines, so 
called Broad Absorption Line Quasars (BAL QSOs). 

Pinally, the separate class of AGN are so called BL Lacerte (Lacs) ob- 
jects. BL Lacs are nuclei of elliptical galaxies that emit highly variable and 
polarized radiation, and have a non thermal optical continuum and a strong 
radio emission. With this group of galaxies, we often connect another type of 
AGN, Optically Violently Variable (OVVs) quasars. The OVVs, in contrast 
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to BL Lacs which have featureless continuum, emit broad optical emission 
lines characteristic for quasars. Both these groups of objects are often re- 
ferred to as one class called blazars. 



2.3 Unified model 



As it was mentioned above, all types of AGN have some common properties. 
It is widely accepted that central engine of an AGN is accretion onto super- 
massive black hole. First observational evidence in favor of a unified model 
of AGN was spectropolarometric observations of Seyfer t 2 galaxy NGC 1068 



(lAntonucci fc Miller! . 119851 ). lAntonucci fc Milleii (119851 ) found broad Balmer 
lines and Fe II emission, characteristic of a Seyfert 1 spectrum, in the polar- 
ized spectrum of NGC 1068. 



In the unified model of AGN (lAntonuccil . 119931 ) both types of Seyfert 
galaxies are intrinsically the same. The difference in spectral characteristics 
appears due to different angle of the visibility of the central regions (see Fig. 

A geometrically and optically thick dusty molecular torus-like structure 
surrounds the central source, as well as the BLR (see Fig. 1). Therefore the 
visibility of the nuclear engine depends on viewing angle with respect to the 
torus. The broad-line polarization in Seyfert 2-s is probably due to electron 
scattering. This picture is consistent with the differenc es observed betwee n 
continuum fluxes of the Seyfert 1 and Seyfert 2 galaxies ( Clavel et al.l . I2OO0I ). 

Following the unified model, the line-of-sights to Seyfert 2 galaxies (or 
Narrow Line Objects) are obstructed by optically thick material correspond - 



ing to hydrogen column densities of Njj > 10 cm (iRisaliti et al.l . Il999l ) 



For column densities < 10 cm~ , photons above a few keV can pene- 
trate the torus creating an un-obstructed view of the nuclear source. In this 
case the source is called Compton thin. For column densities 10^^ > A''^^ > 
10^'* cm~^, only high e nergy X-ray emissio n (tens of keV) can pass through 



the obscuring material (ITurner et al.l . 119971 ). But in the case Nh > 10^^ cm 



even high energy X-rays, above a few tens of keV, are Compton scattered 
and the nuclear source is completely hidden from our direct view. Therefore, 
by observing near the edge-on torus, one can detect narrow line objects, with 
weak (or without) X-ray emission; while by observing near the face-on torus, 
one can detect Blazars or OVVs. The broad line objects can be detected if 
one observes an inclined torus. 
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Figure 1: The scheme of unified model of AGN (Figure courtesy: Vesna 
Borka Jovanovic). 



3 Space-Time Geometry in Vicinity of Su- 
permassive Black Holes 

As it was shown in the previous section, according to unified model of AGN, 
in their heart there is a supermassive black hole with mass up to 10^ Mq, sur- 
rounded by an accretion disk. A black hole is a region of space-time around 
some collapsed mass which gravitational field became so powerful that noth- 
ing (including electromagnetic radiation) could escape fro m its attraction , 



after crossing the certain boundary called the event horizon (IHawkingl . Il988l ) . 

All black holes in nature are commonly classified according to their mass 
as: supermassive black holes (with masses Mbh ~ 10^-10^ Mq), intermediate- 
mass black holes {Mbh ~ 10^ — 10^ ^0), stellar-mass black holes {Mbh < 
10^ Mq), mini and micro black holes {Mbh ^ ^0)- In this chapter our focus 
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will be on supermassive black holes located in the centers of AGN, because 
they are surrounded by accretion disks which X-ray radiation is the main 
subject of this discussion. Space-time geometry in vicinity of a supermassive 
black hole depends on its type (more precisely, on its angular momentum), 
being either non-rotating, where the local space-time geometry is described 
by Schwarzschild metric, or rotating, where the Kerr metric determines the 
geometry of local space-time. 



3.1 Historical background 

The term black hole was introduced in 1969 by the American scientist John 
Wheeler, but the predecessor of the modern idea of black holes emerged in 
the late eighteenth century. It was independently developed by the British 
geologist John Michell in 1783 and the French mathematician and astronomer 
Pierre-Simon Laplace in 1796. This idea was based on Newton's theory of 
gravity and particle theory of light, and is commonly referred to as a dark 
star. A dark star is, according to Michell, a sufficiently massive and com- 
pact star which have such a strong gravitational field that the corresponding 
surface escape velocity equals or exceeds the speed of light. Therefore, any 
light emitted from the surface of such star would be dragged back by the 
gravitational field and we w ould not be able to see it, but we could still feel 



its gravitational attraction (IHawkingl . Il988l ). After the wave theory of light 



was founded, the idea about a dark star was forgotten since, according to the 
wave theory, it was not clear that light could be affected by gravity at all. 

Modern theory of black holes was founded in the twentieth century, more 
precisely in 1915, shortly after Einstein introduced the theory of general 
relativity, when the German physicist and astronomer Karl Schwarzschild 
found a particular exact solution to the Einstein field equations, for the 
limited case of a single spherical non-rotating mass. 

The implications of the general relativity to the stella r evolution wer e 



firstly understood by Subrahmanyan Chandrasekhar in 1928 (jHawkingi . Il988l ) . 
He calculated a maximum limit for the mass of a star, above which the star 
would not be able to support itself (by repulsion of electron degeneracy pres- 
sure which arises from Pauli exclusion principle) against its own gravity, after 
it had used up all its fuel. This limit (~ 1.5 Mq) is nowadays known as the 
Chandrasekhar limit. If a star's mass is less than this limit, it can eventuall y 
stop contracting at a final state known as a "white dwarf flHawkind . [l988l l 



A similar discovery was made in 1932 by the Russian scientist Lev Davi- 
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dovich Landau. He found that there is another hmiting mass (^1—2 Mq), 
for which the corresponding final state is much smaller than a white dwarf 
and which would be supported against its gravity by the exclusion principle 
repulsion between neutrons and protons , rather than be tween electrons. Such 
stars are therefore called neutron stars (iHawkind . Il988l ) . All these results led 
to a surprising conclusion that a star which mass is above the Chandrasekhar 
limit, when it comes to the end of its fuel, could experience a catastrophic 
gravitational collapse to a point with an infinite density. American scientist 
Robert Oppenheimer was the first who explained in 1939 that, according to 
general relativity, stars above ~ 3 Mq would at the end of their lives collap se 
into black holes, due to reasons given by Chandrasekhar (IHawkingl . Il988l ). 

With the discovery of QSOs in the early 60s, the speculations about how 
their enormous luminosities could be produced within such a small region 
occurred. Thermonuclear reactions were not enough efficient, and therefore, 
they were quickly eliminated. Hoyle and Fowler suggested in 1963 that such 
luminosities could be produced during the collaps es and disintegrations o f 
stelar-type objects with masses of 10^ — 10^ Mq (jFerrarese fc Fordl . l2005l ). 
Soon, it became clear that the energy source was gravitational. For example, 
Wheeler suggested a mechanism in which gravitational sin gularity at the cen- 
ter o f a galaxy converted the falling matter into energy (jFerrarese fc Fordl . 
20051 ). Zel'dovich, Novikov and Salpeter proposed in 1964 the release of 



energy through accretion due the growth of a massive object at the cen- 
ter of a galaxy, while Lynden-Bell in 1969 made an attempt to explain the 
phenomena observed in QSOs and Seyfert galaxies in terms of a black hole 
formalism (jFerrarese fc Fordl . |2005| ). A crucial event for the acceptance of 
black holes was the discovery of pulsars by Jocelyn Bell in 1967, because 
it was the clear evidence of the existence of neutron stars, and therefore, 
confirmation of Chandrasekhar limit. The first detection of a solar mass 
black hole came in 1972, when the mass of the rapidly variable X-ray source 
Cygnus X-1 was proven t o be above the maximum allowed for a neutron star 
(jFerrarese fc Fordl . l2005l ) . The first convincing evidence for the existence of 
supermassive black holes is the discovery of t he Fe Kg spectral li ne in X-ray 



spectrum of Seyfert galaxy MGC-6-30-15 bv iTanaka etaP (ll995h FI 



Nowadays, the black hole theory is very well established, mainly due to 
the contributions of numerous authors during the second half of the twentieth 
century, like Roy Kerr, Werner Israel, John Wheeler, Brandon Carter, Roger 



■^See §5.2 for more details 
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Penrose, Stephen Hawking and many others0 

3.2 Schwarzschild metric - non-rotating black hole 

Schwarzschild metric describes the space-time geometry in spherically sym- 
metric gravitational field around a time-steady non-rotating black hole in 
vacuum. The solution of the Einstein field equa tions under these con di- 
tions was found by Karl Schwarzschild in 1916 fISchwarzschildl . Il916al lbl). 



If we denote the Schwarzschild radius which corresponds to a mass M by 

2GM M , , , ^ , , , , 

Rs = — ^ — ~ 3 km (where G is the gravitational constant and c is the 

Mq 

speed of light), then the squared element of space-time interval ds"^ in the 
case of Schwarzschild metric is given by: 



( \ 

1 



dr'-[\-?^\ c^de + r\de'' + sin^ edcf"), (1) 



where t is the time coordinate measured by a stationary clock at infinity 
(which must be distinguished from the proper time r, measured by a clock 
moving with the particle), r is the radial coordinate, 9 is the colatitude (angle 
from North) in radians and is the longitude in radians. 

It is obvious from Eq. ([1]) that Schwarzschild metric depends only on 
the mass M through the Schwarzschild radius Rs- The Schwarzschild metric 
has an infinite space-time curvature (i.e. a singularity) at r = 0. In the 
case when r — >■ cxd this metric is reduced to the Minkowski metric of fiat 
space-time: ds^ = dr"^ — c^dt^ + r'^{d6'^ + sin^ 9d(j?). 

If the complete mass M collapses below its Schwarzschild radius Rs, then 
it will form a spherically symmetric Schwarzschild black hole. A Schwarzschild 
black hole has no charge or angular momentum and it can be distinguished 
from any other Schwarzschild black hole only by its mass. 

To a first approximation, the matter surrounding a black hole may be as- 



sumed to rotate in circular Keplerian orbits with velocity (e.g. lShakura fc Sunyaevl . 
GM 



19731 ): = — — . Therefore, taking into account that a maximum allowed 
R 

rotation velocity is speed of light {v = c), one can obtain so called gravita- 



"^For more detailed historical overview on this topic, see e.g. 



Hawkingl (|l988l ) 
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tional radius Rg: 

GM Rs M ^ 
Ra = —1^ = — ^ 1.5 km, (2) 

which is usually used as a unit for distance around a black hole. 
3.3 Kerr metric - rotating black hole 

Roy Kerr, a New Zealand mathematician, found an exact solution of the 
Einstein field equations of general relativity in the case of gravitational field 
outside an uncharged rotating black hole. The properties of Kerr solution, 
known as Kerr metric, are given i n numerous papers and books (see e.g. 



Chandrasekharl . Il983l : ICarterl . 119731 ) and therefore, here we will present only 
its definition and some basic features, necessary for the further discussion. 

Contrary to the Schwarzschild case, the Kerr metric is no longer spher- 
ically symmetric. In the case of a rotating black hole with mass M and 



angu lar momentum J = aM, this metric is given by (see e.g. iFanton et al. 



19971 ): 



ds^ = -(l- ^) (le-^^ sin^ edtd(t)+^ sin^ edcf'+^dr^+T.de^ (3) 

\ Zj / Zj Zj Za 

where 

E = + cos^ ^, ^ = r'^ + -2Mr, A = (r^ + a^)^ - a^A sin^ 0. (4) 

Above definition is given in Boyer-Lindquist coordinates for c = G = 1. 
As it can be seen from expressions ([3]) and (jl]), in addition to black hole 
mass M, the Kerr metric also depends on its specific angular momentum a. 
Here a < M because any additional angular mo mentum would increase the 



energy of the black hole and, therefore, its mass (IKrolikl . Il999l ). In the case 
of a non-rotating black hole, i.e. for a — > 0, it reduces to the Schwarzschild 
metric. 

By solving the equation A = and taking only a root with + sign, we 
obtain the radius of event horizon of a black hole: Vh = M + V — a?. 
Thus, when a — >• then 2M = Rs (because we assumed that c = G = 
1), whereas when a ^ M (i.e. for a maximally rotating black hole) then 
Th M. 

The minimum allowed radius r^s of a stable circular equatorial orbit, 
or so called marginally stable orbit, is given by the roots of the equation 
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flFanton et all . Il997h : - QMr T Sa^Mr - Sa^ = 0, where the upper sign 
refers to co-rotating orbits, while the lower one to counter-rotating orbits. In 
the case of Schwarzschild metric (a = ) rm.., = 6M, wherea s for a maximally 



rotating black hole (a = M) Vms = M (IFanton et al.l . Il997l ). 



Rotating black holes are surrounded by a region called the ergosphere, 
in which there can be no static observers and the negative energy orbits are 
possible ( iKrolikl . Il999l ). Also, there is a general relativistic effect known as 
frame-dragging, according to which a rotating black hole drags the space-time 
around itself, but such effects will not be discussed here in more details. 



Accretion Disk Around a Supermassive 
Black Hole 



Some observed quasars have luminosities of up to L^oi ~ 10^^ erg/s (ISchneiderl . 
20061 ). The corresponding total energy emitted during the lifetime of such 
quasars can be estimated to E > 3 x 10^^ erg, assuming that their minimum 
age is about 10^ years and that th eir luminosities do not change significantly 
over the lifetime (ISchneiderl . |2006| ) . Only two known mechanisms could pro- 
duce such enormous amount of energy: nuclear fusion and accretion. The 
maximum efficiency e (defined as the mass fraction of fuel that is converted 
into energy, acco rding to E = emc^) of thermonuclear reactions is e < 0.8%. 



Schneiderl (120061 ) showed that this value is too low to explain such high ener- 
gies of > 3 X 10^^ erg. Accretion is the only remaining mechanism which 
can yield larger e. Its maximum efficiency is about e ~ 6% for a non-rotatin g 
black hole and e ~ 29% for a maximally rotating one (ISchneiderl . l2006l ). 
Therefore, AGN release their enormous energy by matter accreting towards, 
and falling into, a central supermassive black hole. 

If the cold matter, which was initially at rest and without magnetic 
field, was subjected to free radial infall it would accret e to the central black 
hole without any energy release or observational effects (IShakura fc Sunyaevl . 
19731 ). However, in the case of AGN the accreting matter has a significant 
angular momentum which does not allow its free infall. At the marginally 
stable orbit of central black hole with mass of lO^M©, the specific angular 
momentum is ~ 1 x 10^^ cm^/s, which is much less even in comparison to a 
typical galaxy wh ere the s peciff ic angular momentum of orbiting material is 
6 X 10^^ cm^/s (iKrolikl . Il999l ). It means that approach of accreting mate- 
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rial toward a black hole requires a loss of the greatest fraction of its initial 
angular momentum. Mechanisms which can contribute to such loss of angu- 
lar moment um are viscosi ty, nonaxisymmetric gravitational forces, magnetic 
forces, etc. f lKrolikl . [T999h . 

Since the orbit of minimum energy for fixed angular momentum in any 
spherically symmetric potential is a circle, the infall of accreting material due 
to loss of its angular momentum will be in form of successively smaller and 
smaller concentric circles. Matter traveling along orbits inclined to each other 
will eventually collide in the plane of intersection, and as a result, the angular 
momenta of different gas steams will be mixed and equalized. Consequently, 
all accreting matter will orbit in a single plane and will have the same specific 
angular momentum at any given radiu s, meaning th at accretion is most likely 
performed through an accretion disk (IKrolikl . 119991 ). 

The accreting material is able to approach the marginally stable orbit 
of central black hole only if there is an efficient mechanism for transporting 
angular momentum outward. Magnetic field that exists in matter fiowing 
into the disk, as well as its turbulent motions, enable angular momentum 
to be transferred outward. Hence, the accretion disk is one of the best can- 
didates for a such mechanism, since the particles in it lose their angular 
momentum due to friction between its adjacent layers and spiral towards 
black hole, releasing their gravitational energy. Part of this energy increases 
the kinetic energy of rotation and th e rest is turned into the the rmal energy 
and irradiated from the disk surface (IShakura fc Sunyaevl . Il973l ) . 

The efficiency of angular momentum transport in the disk is characterized 
by viscosity parameter a < 1, which can be assumed as constant in most 
cases, since the spectrum of disk radiation and its surface temperature do 
not str ongly depend on it, except in the case of a supercritical accretion 
regime (IShakura fc Sunyaevl . Il973l ) . 



4.1 Accretion rate and luminosity of AGN 

The total energy release and the spectrum of emitted radiation are mainly 
determined by the rate of matter inflow into the disk on its outer boundary, 
i.e. by its accretion rate which is usually denoted by M. If the accretion 
converts matter to radiation with flxed radiative efficiency 1] (in rest-mass 
units), then a characteristic scale for the accretion rate is the Eddington 
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accretion rate (jKrolikl . Il999l : IShakura fc Sunyaevl . Il973l ): 

0.06 M Mq 



3 X 10 



-8 



V Mq yr 



(5) 



for witch the total release of energy in the disk L 
Eddington luminosity: 

,38 M erg 



1.51 X lO'^ 



TjMc^ is equal to the 



(6) 







Eddington luminosity defines a critical luminosity for any given mass, beyond 
which the radiation force overpowers gravity. 

Observed AGN have luminosities from ~ 10^^ to ~ 10^^ erg/s, which 
means that t heir central b lack holes must have masses from 10^ to 10^ M©, 
respectively (Krolik, 19991 ). If we assume that accretion in AGN occurs with 
rate of M ~ 0.5 Mq/yt and with radiative efficiency r] ~ 0.1, the resulting 
luminosity would be ~ 3 x 10^^ erg/s, which i s the raiddle of the observed 
luminosity distribution of low redshifted AGN (jKrolikl . Il999l ) . 

Actual accretion rates M of matter inflow into the accretion disks of AGN 
could be many times less or higher than the critical value of the Eddington 
accretion rate. At essentially subcritical accretion rates (M -C Me), the 
maximal surface temperatures Tg are on the order of 10^ — 10^ K in the inner 
regions of the disk, fro m which the most of energy is released, mainly in UV 
and soft X-ray bands (IShakura fc Sunyaevl . Il973l ). When the accretion rate 
increases, the luminosity also raises, as well as the effective radiation tem- 
perature. At accretion rates comparable to Me, the accretion disk becomes 
a powerful source of X-ray radiation with the effective radiation tempera- 
ture Tr from 10'' to 10^ K. In a strongly supercritical regime of accretion 
(M 3> Me), the luminosity becomes fixed at the Eddington critical limit Le 
and the most o f energy is irradiated frorn accretion disk in UV and optical 
spectral bands (IShakura fc Sunyaevl . Il973l ) . 



4.2 Standard model and spectral distribution 

By the term "standard Newto nian model of accret i on di sk", here we will 
assume the disk model given by lShakura fc Sunyaevl (Il973l ). This model was 
originally developed to describe the accretion disks around stellar sized black 
holes in the binary systems, but with certain modifications it could be also 
applied on accretion disks around supermassive black holes in AGN. 
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The standard Newtonian model of accretion disk is based on the supposi- 
tion that the released gravitational energy is emitted as a multitemperature 
blackbody radiatiori , wher e the surface temperature profile is given by (see 
Shakura fc Sunyaevl . Il973l ): 



TJr) 



3 GM 



57r err- 



M 1 - 




1/4 



(7) 



where a is the Stefan constant, G is the gravitation constant, M is the mass 
of the central black hole, M is the accretion rate and rj„ is the inner radius 
of the accretion disk. In above Equation and t hroughout this section we will 
use the dimensionless disk radius r defined as (jPopovic et all l2006al ) : 



R 



1 m 



Mq r 

6 Ra ~ 6 GM ~ M 9 km' 



(8) 



where R is disk radius expressed in gravitational radii Rg. 

In the standard model of accretion disk, accretion occurs via an opti- 
cally thick and geometrically thin disk. The effective optical depth in the 
disk is very high and p hotons are close to thermal equilibrium with electrons 



(iPopovic et al.l . I2006al ). The spectrum of thermal radiation emitted from ac- 



cretion disk surface depends on its structure and temperature (and therefore 
on the distance to the black hole), and can have several distributions such 
as ( Shakura &: Sunyaev . 19731 ): 



1. Planck distribution in the outer regions of the disk, where free- free 
and free-bound processes, as well as absorption in the lines of heavy 
elements broadened by the gas pressure, give the main contribution to 
the shape of emitted spectrum: 



Fix) 



27Th fkT 



X 



- 1' 



(9) 



hu 



where x = - — , h is Planck and k is Boltzman constant, 

more complex spectrum which passed through the homogeneous or ex- 
ponentially varying medium in intermediate regions of the disk where 
Thomson scattering dominates and 
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3. Wien distribution in the inner regions where the comptonization pro- 
cesses strongly affect the shape of emitted spectrum: 



IX 



(10) 



The surface temperature has different radial distributions in different 
parts of accretion disk and results in the multicolor black body spectrum. 
In the inner regions of accretion disks of AGN where their X-ray radiation 
is gen erated, the rad i al distr ibution of surface temperature has the following 
form (jPopovic et al.l . l2006bl ) : 



To r-3/2(i_ ^-1/2)4/5 



(11) 



where the temperature constant Tq is chosen so that the corresponding ef- 
fective radiation temperature is Tr = 10^ — 10^ K. This modification of stan- 
dard disk model is made, because the surface temperature given by Eq. ([7]) 
could not be successfully applied in the case of accretion disks of AGN. The 
distribution of the temperature as a function of the accretion disk radius, ac- 
cording to Eq. ( fTTj) , is presented in Fig. [2] for two differen t values of angular 
momentum a of central black hole (jPopovic et al.l . l2006al ) . 



a.0.002 
a.0.998 




Figure 2: The distribution of the temperature as a function of the disk radius 
R (in the direction normal to the rotation axis, as seen by distant observer), 
given for two different values of angular momentum a of central black hole 
(jPopovic et al.l . l2006al ). Negative values of R correspond to the approaching 
and positive to the receding side of the disk. 



One can see from Fig. [2] that the temperature depends not only on disk 
radius, but also on its other parameters, such as e.g. angular momentum a. 
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Emitters located at different radii in accretion disk have different tempera- 
tures and therefore make different contributions to the observed ffux. Also, 
it is noticeable that the temperatures at smaller radii of approaching side of 
the disk are significantly higher than those on its receding side. 

Intensity of radiation emitted from the total disk surface can be obtained 
by integration of all local spectra emitted from the above regions: 

Iu = ^n j F,[Ts{r)]rdr, (12) 

''in 

where rj„ and Tout are inner and outer radii of accretion disk, respectively. 







4. 





Figure 3: Magnetohydrodynamic simulations of an accretion disk as seen by a 
distant observer at the following inclination angles (lArmitage and Reynolds! . 
2003f ): 5° (upper left), 30° (upper right), 55° (lower left) and 80° (lower 
right). The disk rotation is in counterclockwise direction. 



Magnetohydrodynamic simulations of an accretion disk as se en by a dis- 

tant observer at four different inclinations are presented in Fig. [3] (lArmitage and Reynolds 
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20031 ). The most distinctive details in this figure are especially bright regions 



in form of arcs within the turbulent flow, which trace out the photon tra- 
jectories close to the radius of marginally stabile orbit. These arcs are much 
brighter on approaching side of the disk due to higher surface temperature, 
as shown in Fig. [2], and also due to Doppler boosting and relativistic beaming 
which both enhance the flux observed in the direction of motion and diminish 
the flux in the opposite directionjf] 



4.3 Structure and emission 

An accretion disk around a supermassive black holes in the center of AGN 
extends from the radius of marginally stable orbit Rms to the several thou- 
sands of gravitational radii. According to its radia tion emitted in differe nt 



spectral bands, it can be stratifled in several parts (I Jovanovic et al.l . l2008l ): 



innermost part close to the central black hole which emits X-rays and 
which extends from the radius of margi nally stable orbit -Rm., to th e 
several tens of gravitational radii (see e.g. iBallantyne fc Fabianl . l2005l ). 



2. central part ranging from 
radiation and 



100 Rg to ~ 



1000 Eg which emits UV 



3. outer part extending from several hundre ds to several thousands R 



from which the optical emission is coming (lEracleous fc Halpernl . Il994 



20031 ). 



Emissivity of the disk has an important role in the line and continuum 
shapes. In the rest frame of the emitting material (i.e. in disk frame), 
emissivity e( r) is deflned as e n ergy emitted per unit proper time per unit 



proper area (jPabrowski et al.l . 119971 ). In the same fram e, it is related to 



the e mitted intensity /(r, i^e) and emitted flux F{r, Vg) by (jPabrowski et al 
19971 ): 



Fir, 



Trl(r, Vf. 



e{r)8{vg - z/q) 



(13) 



under assumption that the line emitted at frequency z/q can be approximated 
by a 5 function. 

Although the standard model of the accretion disk does not predict the 
power-law for the disk emissivity, such law is usually accepted in the case 



^More details about the last two relativistic effects can be found in e.g. iKrolikI (Il999h 
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of the hard X-ray emission from t he inner pa rts of AGN accretion disks 
(jFabian et al.l . Il989l : iNandra et all 119971 . Il999l ). In that case the surface 
en aissivity is assumed to vary with radius a s power law with emissivity index 
p fiFabian et al.l . Il989l : IPanton et all Il997h : 

^0 



eir 



(14) 



where eo is the emissiyity constant. iFabian et al.l (1l989l ) found that emissiyity 
index p varies between and 3 for the Fe Ka line emitting region of Cygnus X- 
1, where the most probable value is p ~ 2. Some authors found that models 
describing the disk emissivity as a broken power law, rather than a single 
power law, achi eve significantly better statistic al fits of the X-ray emission in 
the Fe Ka line ( iBrenneman &: Reynolds! . |2006| ) . One such broken power law 
between s ome inner radius rjn and outer radius Vout with break radius rtr is 
given by (IBrenneman fc Reynolds! . |2006| ) : 

0, 



eir 



{r/rbr 



-a2 



0, 



r < rin 
^in r <^ r}y^ 

r > rout 



(15) 



where ai and a2 are some emissivity indices. Using this law IBrenneman fc Reynolds 
fl2nnfi! ) iitted the Fe Ka line of MCG-6-30-15, assuming a maximally spin- 
ning (a = 0.998) black hole, and found that the line emissivity followed such 
dependence for ai = 4.5 — 6 and a2 ~ 2.5, where the break radius was 
rbr ~ 6 Rg. 

The X-ray continuum can be also modeled e.g. by the following time- 
independent intri nsic emissivity o f the disk as a funct i on of i ts radius r and 
photon energy E (IDovciak et al.l . |2004J : iPopovic et al.! . l2006a! ) : 

/(£;,r) ~ E"^ X r-°, (16) 

where F is photon index and a is radial emissivity power law index of a 
continuum emitting region. According to some studies of observed X-ray 
spectra of AGN, the values of F and a are estimated to be about 1.5 and 
2.5, respectively f Povciak et al. . 20041 : Popovic et al. . 2006a ). 

In the case of the disk outer parts (i.e. for optical emitting region), the 
black-body emissivity law is assumed and, according to the Eg. ([9]), emitted 
intensity is given by Planck function (IPopovic et al.l . l2006a! ): 

2E^ 1 
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However, in the inner parts of the accretion disk Planck function cannot 
be used properly. In these regions so called "modified" black-body emissiv- 
ity law is applied. According to the Eg. (fTOl) . the corresponding emitted 
intensity is given by (jPopovic et al.l . l2006al ) : 

i{e,t; ' - ^ 



oc x^e 



X 



a 
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Figure 4: Shapes of the X-ray continuum in 1 - 10 keV energy band, emitted 
from an accretion disk with the inner and outer radii equal to Rin = Rms 
and Rnut = 80 Rg, arou nd a supermassive black hole with mass M = 10^ Mq 



(iPopovic et al.l . l2006al ). The continuum shapes are obtained assuming the 



following three emissivity laws (all normalized to their maximum values): 
the power law defined by Eq. ([T3j) with emissivity index p = 2.5, the black- 
body and "modified" black-body emissivity laws defined by eqs. f|T7|) and 
(|T8l) . respectively, assuming the radial distribution of the surface temperature 
given by Eq. ffTTj) . 



A comparison between the shapes of the X-ray continuum (1-10 keV) 
emitted from an accretion disk around a supermassive black ho le, obtained as- 



sumin g three different emissivity laws, is presented in Fig. H] (IPopovic et al 



2006al ). From Fig. [Hit is clear that X-ray continuum shape strongly depends 
on adopted emissivity law. 
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5 Supermassive Black Holes and X-ray Emis- 



sion 

A disk geometry for the X-ray emitting regions of AGN could be assumed, 
since the unifie d model of AGN inc ludes a supermassive black hole fed by an 
accretion disk. iFabian et al.l (119891 ) calculated spectral line profiles for radi- 
ation emitted from inner parts of a ccretion disks and la ter on such features 
of Fe Ka lines were discovered by iTanaka et al.l (119951 ) in Japanese ASCA 
satellite data for Seyfert galaxy MGC-6-30-15. The assumption of a disk ge- 
ometry for the distribution of the X-ray emitters in the centra l parts of AGN 



is also supported by the spectral shape of the Fe Ka line ( iNandra et al. 



19971 . 11999I . I2OO7I ). Moreover, a bump in the UV spectra of AGN is present 
very often, indicat ing that the UV and optical continua also originate from 
an accretion disk (jjovanovic et al.l . 120081 ). 

We should note here that probably most of the X-ray emission in the 1- 
10 keV energy range originates from inverse Compt on scattering of photon s 



from the disk by electrons in a tenuous hot corona (jJovanovic et al.l . |2008| ) . 



Proposed geometries of the hot corona of AGN include a spherical corona 
sandwiching the disk and a patchy corona r aade of a few compact regions 
covering a small fraction of the disk (iMalzad . 120071 1. On the other hand, it is 
known that part of the accretion disk that emits in the 1-10 keV rest-frame 
band (e.g. the region that emits the continuum Compton reflection compo- 
nent and the fluorescent emission lines) is very compact and may contribute 
to the X-ray variability in this energy range. 

A typical X-ray sp ectrum of AG N (see Fig. [5]) is composed from the 
following components (iFabianl . |2006| ) : 



- an underlying power-law continuum due to thermally Comptonized soft 
photons, 

- a soft excess at low energies below 1 keV due to thermal (black-body) 
emission from an optically-thick accretion disk, 

- a fluorescent/recombination Fe Ka line and 

- a Compton hump due to X-ray reflection from the disk. 



In the following text we will discuss in more details power law component 
of X-ray continuum and the Fe Ka spectral line. 
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5.1 X-ray continuum of AGN 



AGN are powerful sources of X-ray radiation in the continuum from 0.1 to 
100 keV, which contains two components (see Fig. [6]): the soft X-rays with a 



(Fabian et aL. 


1989; 


Fabian . 


2006) 



Therefore, the observed continuum flux 
is very often fitted with one or two black-body components in the soft X-rays, 
in addition to a power law component in the hard X-rays. It is believed that 
both, soft and hard components mostly arise from the inner region of AGN 
disk, close to the central supermassive black hole (~ 10 Rs)- The first com- 
ponent is probably formed in the accretion disk, which is also a strong source 
of the soft UV and optical photons, or at least in cold (T < 10^ K) accreting 
gas clouds. The second one is caused by high-energy (most likely relativistic) 
electrons in hot corona above the disk, when they multiply inverse- Compton 
scatter some of the low-eriergy UV an d optical photons from the disk to X- 



ray energies (IFabian et al.l . Il989l . l2000l ). The resulting hard X-ray power law 



component irradiates the accretion disk and produces a reflection component 
which causes the observed spectrum to flatten a bove 10 keV (see Fi g. [6]), as 



Compton recoil reduces the backscattered flux (IFabian et al.l . |2000| ). Apart 



from its own thermal radiation, the cold gas around a black hole is also ir- 
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radiated by these hard X-ray photons, producing different spectral features 
through photoelectric absorption, fluor escence (responsible for occurrence of 



fluorescent iron lines) and scattering (IFabian et al.l . Il989l ). Therefore, the 



observed X-ray spectra of AGN consist of an ir on line spectrum superposed 



on a power law continuum (IFabian et al.l . |2000| ). 



6- 




1 2 
Energy (keV) 



Figure 6: The two-component model of the X-ray continuum in the case of 
Seyfert 1 galaxy IH 0707-49 5, composed from a power law (PLC) and a 
reflection (RDC) component (IFabianl . |2006| ). 



Fluctuations of the X-ray radiation on timescales from s everal parts of 



19991). 



li p 

an hour until several days are a common property of all AGN (iKrolikl . 
Such fast variations confirm assumption that X-ray radiation is emitted from 
a v ery cornpact r egion in the center of AGN. 

FabianI (120061 ) showed that a two-component model, consisting of a power- 



law and a reflection component, could better describe the observed variability 
of X-ray spectra in the case of several bright Seyferts than a simple power-law 
model. This two-component model is presented in Fig. [61 



5.2 Fe Ka spectral line 



Contrary to the UV and optical spectra of AGN, there are very few strong 
spectral lines in their X-ray spectra. The most important one is the Fe Ka 
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line. Iron abundance in accretion disks of AGN is sufficient to produce very 
strong emission in this line. The CCD detectors on Japanese ASCA satellite 
were the ffist instruments with sufficient spectr al resolution and sensitivity in 
the X-ray band, by which iTanaka et al.l (119951 ) obtained the first convincing 
proof for the existence of the Fe Ka line in AGN spectra. This discovery was 
made after four-day observations of Seyfert 1 galaxy MCG-6-30-15 (see Fig. 

ED- 




Energy (keV) 



Figure 7: The profile of the Fe K a line from Seyfert 1 galaxy MCG-6-30-15 
observed by ASCA SIS detector (ITanaka et al.l . 119951 ) and the best ffi (blue 
solid line) obtained by a model of the accretion d isk in Schwarzschild metric, 
extending between 3 and 10 Schwarzschild radii (IFabian et al.l . 119891 ). Image 
Credit: Tanaka et al.. Copyright: Nature, 375, 659, 1995. Image generated 
by Dr. Paul Nandra NASA/GSFC. 

The fluorescent/recombination iron Ka line is an important indicator of 
accreting flows around compact objects, because it is produced in inner parts 
of their accretion disks. At the same time, it is the strongest line of the X-ray 
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radiation, and it can be found in the spectra of all types of accreting sources: 
binary black hole and neutron star systems, cataclysmic variable stars and 
AGN. 

The Fe Ka line is produced when plasma is subjected to the influence of 
the hard X-ray radiation so that one of the two i^'-shell [n = 1, where n is 
the principal quantum number) electrons of an iro n atom (or ion) is e jected 



following the photoelectric absorption of an X-ray fiFabian et all. l2000[i_^_Th e 



2000h . 



threshold for the absorption by neutral iron is 7.1 keV (iFabian et al.l . 
The resulting excited state decays when an L-shell (n = 2) electron drops 
into the ii'-shell, releasing 6.4 keV of energy. This energy is either emitted as 
an emission-line photon (34% probability) or internally absorbed by another 
electron (66% probability) which is consequently ejected from the iron ion 
(Auger effect). 

The fluorescent yield (i.e. the probability that photoelectric absorption 
is followed by fluorescent line emission rather than the Auger effect) is a 
weak functiori of the ionization state from neutral iron (Fe I) up to Fe XXIII 
( Fabian et al.l . 2000 ). For lithium-like iron (Fe XXIV) through to hydrogen- 
like iron (Fe XXVI), the lack of at least two electrons in the L-shell means that 
the Auger effect cannot occur. For He and H-like iron ions, the line is pro- 
duced by the capture of free electrons (recombination) and the equivalent flu- 



oresc ent yield is high and it depends on the plasma conditions ( iFabian et al. 
2000f ). 

For the neutral iron, the Fe Ka li ne energy is 6.4 ke V (more precisely 



there are two components of the line ( Fabian et al. . 2000l ): Fe Kai at 6.404 



and Fe Ka2 at 6.391 keV), while in the case of ionization, the energy of both 
the photoelectric threshold and the Fe Ka line are slightly increased. Even 
for such high ionization states of He and H-like iron ioi is, the Fe Ka line 
energy is increased only to 6.7 and 6.9 keV, respectively (iKroliki . Il999l ). 

Fe Ka line is pretty narrow in itself, but in case when it originates from 
a relativistically rotating accretion disk of AGN it becomes wider due to 
kinematical effects, and also its shape (or profile) is changed due to Doppler 
boosting and gravitational redshift. Such broadening of the line is very of- 
ten observed in spectra of Seyfert galaxies and is one of the main evidences 
for the existence of a relativistic accretion di sk wh i ch ext ends deeply in the 
gravitational field of the central black hole ( Zyckil . 2004). If the line origi- 
nated from an arbitrary radius of a nonrelativistic (Keplerian) accretion disk 
it would have a symmetrical profile (due to Doppler effect) with two peaks: 
a "blue" one which is produced by emitting material from the approaching 
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side of the disk in respect to an observer, and a "red" one which corresponds 
to emitting material from the receding side of the disk. The widest parts of 
the Fe Ka line arise from the innermost regions of the disk, where the rota- 
tion of emitting ma terial is the fastest. Using ASCA satellite observations, 
Nandra et al. fll997h found that, in case of 14 Seyfert 1 galaxies. Full- Widths 



at Half-Maximum (FWHM) of their Fe Ka lines correspond to velocities of 
~ 50,000 km/s. In some cases (e.g. for Seyfert 1 g alaxy MCG-6-30-15 ), 
FWHM velocity reaches 30% of speed of light (see e.g. iNandra et al.l . 120071 ). 
It means that in the vicinity of the central black hole, orbital velocities of 
the emitting material are relativistic, causing the enhancement of the Fe Ka 
line "blue" peak in regard to its "red" peak (relativistic beaming). Tak- 
ing into account the integral emission in the line over all radii of accretion 
disk , one can obt a in the line with asymmetrical and highly broadened pro- 
file fiFabian et all . l2000f ). The "blue" peak is then very narrow and bright, 
while the "red" one is wider and much fainter (see Fig. [7]). Besides, the 
gravitational redshift causes further deformations of the Fe Ka line profile 
by smearing the "blue" emission into "red" one. Since the observed Fe Ka 
line profiles are strongly affected by such relativistic effects, they represent a 
fundamental tool for investigating the plasma conditions and the space-time 
geometry in the vicinity of the supermassive black holes of AGN. 

One of the important features of the Fe Ka line is variability of both, 
its shape and intensity. Observed variations of this line are surprisingly less 
than those of the hig h energetic c ontinuum, which is assumed to give rise 
to the line emission (IZyckil . I2004J ). Also, it seems that there is a lack of 
corresponding line response to the continuum variations on time scales from 
several minutes t o several day s, or that these line and continuum variations 
are uncorrelated (IZyckil . 120041 ). In §6 of this chapter we will pay attention to 
some possible causes of such behavior of the Fe Ka line and X-ray continuum. 



5.3 Modeling of X-ray emission using ray-tracing in 
Kerr metric 



The disk emission can be analyzed by nu merical simulations, based on so 



called ray-tracing method in Kerr me t ric (|Bao et al.l . 1 19941 : iBromley et al 



19971 : iFanton et al.l . 119971 : ICadez et al.l . Il998l ). taking into account only pho- 
ton trajectories reaching the observer's sky plane. In this method one divides 
the image of the disk on the observer's sky into a number of small elements 
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(pixels). For each pixel, the photon trajectory is traced backward from the 
observer by following the geodesies in a Kerr space-time, until it crosses the 
plane of the disk (see Fig. [8]). Then, the flux density of the radiation emit- 
ted by the disk at that point, as well as the redshift factor of the photon 
are calculated. In that way, one can obtain the color images of the accretion 
disk which a distant observer would see by a high resolution telescope. The 
simulated line profiles can be calculated taking into account the intensities 
and received photon energies of all pixels of the corresponding disk image. 
All illustrations of the accretion disk and the line shape in this chapter are 
obtained using such numerical simulations. Here we will briefly describe a 
pseudo-analytical approach of ray-tracing proposed by lCadez et al.l (Il998l ). 




flat space-time 



Figure 8: Schematic illustration of the ray-tracing method in the Kerr metric, 
showing a light ray emitted from some radius of accretion disk in coordinate 
system defined by a rotating black hole with angular momentum a, and ob- 
served at a pixel with coordinates (impact parameters) a, f3 on the disk image 
in the observer's reference frame (Figure courtesy: Vesna Borka Jovanovic). 



This method is based on the pseudo-analytical integration of the geodesic 
equations which describe the photon trajectories in the general case of a 
rotating black hole having some angular momentum J, which gravitational 
field is therefore described by the Kerr metric (see §3.2 for more details). 
The Kerr metric depends on the angular momentum normalized to the mass 
M of black hole: a = J/Mc, < a < M. 
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A photon trajectory in the Kerr metric can be described by three con- 
stants of motion (the energy at infinity and two constants related to the an- 
gular momentum, respectively) wh ich have the follow ing forms when natural 
units c = G = M = 1 are assumed (ICadez et al.l . ll998l ): E = —pt, A = and 
Q = Pq — + A'^cot'^9 . Here, (r, 9, (f), t) are the usual Boyer-Lindquist 

coordinates and p is the 4-momentum. As the trajectory of a photon is inde- 
pendent on its energy, it may be expressed using the two dimensionless pa- 
rameters X = A/E and q = Q^l"^ jE which are very simply related to the two 
impact parameters a and /3 describing the apparent position on the observer's 

celestial sphere: a = —— — - — and P = ± (g^ -|- cos^ 9obs — cot^ 9obs) ^ , 



sin 9, 



obs 



— 

d9 

In order to find the photon trajectories (null geodesies) which originate 
in the accretion disk at some emission radius and reach the observer at 



obs 



infinity, one must solve the following integral equation (jCadez et all Il998l ): 



± 



dr 



± 



d9 



{9,x,qy 



(19) 



aXy - A[{\-ay + q^]. 
2 ' a^cos^9- \^cot^9. 



(20) 



/2(r,A,g) 

e(^,A,g) = r 

Above integral Equation (fT9|) can be solved in terms of Jacobian elliptic 
functions, and therefore it is a pseudo- analytical in t egrati on. For the exact 
expressions of the solutio ns, see e.g. ICadez et al.l (1l998l ). or more recent 
paper from iLi et al.l (120051 ) . 

Due to relativistic effects, photons emitted at frequency Ugm will reach 
infinity at frequency Uobs, and their ratio determines the shift due to these 
effects: g = The total observed flux at the observed energy Eobs is given 



by fiFanton et al.l . 119971 ): 



Fobs {E, 



obs ) 



e {r)g'^S {Eobs - qEq) dl 



(21) 



image 



where e (r) is the disk emissivity, d'E is the solid angle subtended by the disk 
in the observer's sky and Eq is the rest energy. 

A simulated accretion disk image is obtained in the following way: 
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1. values of the following input parameters are specified: inner and outer 
radii and Rout) of the disk, angular momentum a of the central 
black hole, observer's viewing angle (disk inclination) 9obs (also, often 
denoted by i) and parameters defining the disk emissivity 

2. constants of motion A and q are calculated for each pixel on imaginary 
observer's photographic plate (i.e. for each pair of impact parameters 
a and P) 

3. geodesic Equation (fT9l) is integrated for each pair of A and q 

4. values of shift due to relativistic effects g and observed fiux Fobs are 
calculated 

5. pixels on imaginary observer's photographic plate are colored according 
to the value of shift g and a simulated disk image is obtained. 

The simulated line profiles can be calculated from the corresponding disk 
images by binning the observed flux at all pixels over the bins of shift g. 
The examples of simulated disk images obtained in such way are presented 
in left panels of Figs. [9] - [TTl and the corresponding simulated line profiles 
are presented in the right panels of the same figures. 



5.4 Observational effects of strong gravity in the vicin- 
ity of supermassive black holes 

In general, black holes have three measurable parameters (not including the 
Hawking temperature): charge, mass (and hence gravitational field) and 
angular momentum (or spin). In the case of supermassive black holes of AGN, 
only the latter two are of sufficient importance because they are responsible 
for several effects which can be detected in the observed Fe Ka line shapes 



(jJovanovic and Popovid . l2008af ). 

In order to study the size of the Fe Ka line emitting region, as well as its 
location in the disk, one can assum e that the line is emitted from a region in 



form of a narrow ring. For example, iJovanovic and Popovic fl2008ah assumed 



a line emitting region with width equals to 1 Rg, located between: a) = 6 
Rg and Rout = 7 Rg and b) Rin = 50 Rg and Rout = 51 Rg. These two cases 
are presented in Fig. [91 

From Fig. [9] one can see how the Fe Ka line profile is changing as the 
function of distance from central black hole. When the line emitters are 
located at the lower radii of the disk, i.e. closer to the central black hole, 
they rotate faster and the line is broader and more asymmetric (see Fig. IH] 
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X (in Rg) 




X (in Rg) 



9 = E„t,/E„ 



Figure 9: Left: illustrations of the Fe Ka line emitting region in form of 



narrow ring with width 
(top) and Rin 



IRq, extending from: 



50 to i?, 



V 

mil 



6 Rg to Rout 



7R„ 



51 Rq (bottom). Right: the corresponding 



Fe Ka line profiles (jjovanovic and Popovid . l2008al ) . 



top-right). If the line emission is originating at larger distances from the black 
hole, its emitting material is rotating slower and therefore the line becomes 
narrower and more symmetric (see Fig. [9] bottom-right). In majority of 
AGN, where the broad Fe Ka line is observed, its profile is more similar to 
the modeled profile as obtained under as sumption that the li r ie emitters are 



located close to the central black hole (iTanaka et al.l . Il995l : iNandra et al 
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20071 : IJovanovic and Popovid . l2008al ) . 



Angular momentum or spin of the central supermassive black hole of 
AGN is a property of the space-time metric. To demonstrate how it affects 
the observed line profiles we will now assume that the Fe Ka line emitting 
region extends between the following inner and outer radii: i?^ = Rms and 
Rout = 20 Rg. We will analyze the following two cases for accretion disk 
inclination in both Schwarzschild and Kerr metrics: (i) i = 35° and (ii) 
i = 75°. In the case of Schwarzschild metric we have a stationary black 
hole and hence, angular momentum (normalized to the mass of black hole) 
is a = 0. In Kerr metric (i.e. for a rotating black hole), it can take any value 
from the [0, 1] range, but in these two examples we will assume an almost 
maximally rotating black hole with a = 0.998. 

Illustrations of an accretion disk and the corresponding Fe Ka line shapes 
in the first case are presented in Fig. [TOl As one can see from Fig. [TOl the 
red peak of the Fe Ka line is brighter in case of almost maximally rotating 
black hole, but at the same time it is also more embedded into the blue peak 
wing and therefore less separable from it. 

Fig. [m contains illustrations of the line emitting regions and the corre- 
sponding line shapes in the case of a highly inclined disk {i = 75°). Here, the 
line profiles are broader than in the first case, mostly due to higher inclina- 
tion. As it can be seen from Fig. [Ill iii the case of the Kerr metric, the red 
peak of the line is again more embedded into its blue peak wing (as in the first 
case) and it confirms that this effect can be most likely attributed to angular 
momentum (IJovanovic and Popovid l2008a] ). Consequently, angular momen- 
tum of the central black hole has significant influence on the line shape which 
supports assumption that the line originates from the innermost part of ac- 
cretio n disk, close to the central black hole (see e.g. iBallantyne fc Fabianl . 



20051 ). This fact can be used for estimation of angular momentum of the 



central black hole in observed AGN. 

Above simulations of the strong gravitational field influence on the Fe 
Ka line show that such effects can be detected in the observed line shapes 
and therefore, comparisons between the observed and modeled Fe Ka line 
profiles (see e.g. Fig. [7]) can bring us some essential information about the 
strong gr avitational fie l d in t he vicinity of central supermassive black holes 
of AGN fiNandra et all . l2007f ). 
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X {^n Rg) 



9 = E„b,/Ee. 



Figure 10: Illustrations of accretion disk (left) and the corresponding Fe Ka 
line profiles (right) in the case of Schwarzschild (top) and Kerr metric with 
angular momentum a = 0.998 (bottom). The disk inclination is i = 35° 
and its inner and outer radii ar e Rm = Rms and Rout = 20 R^, respectively 
( Jovanovic and Popovic . 2008al ). 




Figure 11: The same as i n Fig. [T0| but for a highly inchned disk with i = 75 
(jjovanovic and Popovid . l2008al ). 
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6 Variability of X-ray Emission Around a Su- 
permassive Black Hole 

Rapid and irregular variability of the observed X-ray emission in the line, as 
well as in the continuum, is a common property of all AGN. This variability 
could be due to disk instability, reflecting in perturbations of the disk emis- 
sivity, or it could be caused by some external effects, such as gravitational 
microlensing and absorption by X-ray absorbers. We developed a model of 
perturbations of the disk emissivity, a model of absorption region and three 
models of gravitational microlensing. In the following text we will pay more 
attention to all three mentioned causes of the X-ray variability in AGN. 



6.1 Perturbations of disk emissivity 



In some cases the observed X-ray variability of AGN cannot be explained by 
the standard model of accretion disk. For example, in addition to the stable 
6.4 keV core of the Fe Ka line, a variable "red" feature of the line at 6.1 
keV is also detected in X-ray spectrum o f Seyfert galaxy NGC 3516 (see Fig. 



[T2D, observed by XMM-Newton satellite (llwasawa et al.l . |2004| ). This feature 



varies system atically in the flux at i ntervals of 25 ks and in energy between 5.7 
and 6.5 keV. llwasawa et al.l (120041 ) found that the spectral evolution of the 
"red" feature agrees well with hypothesis of an orbiting spot in the accretion 
disk. 

Many processes in the accretion disk may lead to perturbations in its 
emissivity, such as self gr avity, disk-star collisions and baroclinic vorticity 
fiFlohic fc Eracleousl . l2008h . Different models of emissivity perturbing regions 
can be used to describe the observed variability in the Fe Ka line profiles 

and in tensities, like stochastically perturbed one given bv lFlohic fc Eracleoud 

( 20081) . but here we will present a bright spot model given by I Jovanovic and Popovic 
(l2008bl ). In this model a modification of the power-law disk emissivity is pro- 
posed in order to explain the observed profiles. The following emissivity law 
of the disk is assumed: 

/ ( fx - Xr>V (y -y^YW 



£i{xp,yp) = e{r{xp,yp)) 



1 + Ep - e 



(22) 



V / 

where ei{xp,yp) is the modified disk emissivity at the given position {xp,yp) 
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Time (s) 



Figure 12: Light curves of Seyfert galaxy NGC 3516 (llwasawa et al.l . l2004j ) 
for: the 0.3 - 10 keV band (top), the Fe Ka hne red feature (middle) and 
the 6.4 keV line core (bottom). 



of perturbing region (in gravitational radii Rg), e{r{xp,yp)) is the power-law 
disk emissivity at the same position, Sp is emissivity of perturbing region and 
{wx,Wy) are its widths (also in Rg). 

An example of the shape of the perturbed emissivity for an accretion 
disk in Schwarzschild metric, as well as the corresponding perturbed and 
unperturbed Fe Ka line profiles, are presented in the left and right panels of 
Fig. [ini respectively. 

As one can see from Fig. [131 when perturbation moves along receding 
side of the disk (positive direction of a;-axis), the perturbing model affects 
only "red" part of the line flux, while the "blue" one stays nearly constant, as 
well as the line core. Therefore, this bright spot model of perturbing region 
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Figure 13: Left: shape of perturbed emissivity of an accretion disk in 
Schwarzschild metric for the following p arameters of perturbing region : 



20 Rg, Up = and 



7 Rg ( Jovanovic and Popovic . 2008bl ). 



Right: the corresponding perturbed (dashed line) and unperturbed (solid 
line) Fe Ka line profiles. 



could satisfactorily explain the observed variations of the Fe Ka line flux. 

Under assumption that perturbation moves by speed of light c, one can 
also calculate time tp [s] that corresponds t o the current position (xp,yp) o f 
perturbation, using the following expression ( Jovanovic and Popovic . 2008bl ): 



^ r n _ r{xp,yp) [Rg] ^ r^^y ■ GMbh 
^ c[m ■ 

T (x , y ) 

where r^y = ^ , G is Newton's gravitational constant and Mbh is the 

Rg 

mass of central black hole. 

Using the time tp [s] we are now able to obtain simulated light curves, 
produced as perturbation moves along the accretion disk. An example of 
such light curves, corresponding to displacement of perturbation along the 
receding side of the disk, is given in Fig. [TH From this figure it can be seen 
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Figure 14: The simulated light curves in case of perturbing region 
displacement along the recedi ng side of the disk and without rotation 
( Jovanovic and Popovic . 2008bl ). Light curves correspond to the following 
spectral bands: total flux (black) to 0.1 - 12.8 keV, "red" to 0.1 - 6.1 keV, 
"core" to 6.1 - 6.7 keV and " blue" to 6.7 - 12.8 keV. 



that displacement of perturbing region results in variations of only "red" light 
curve (0.1 - 6.1 keV), while the "blue" one (6.7 - 12.8 keV) and the line core 
(6.1 - 6.7 keV) stay nearly constantly Thes e variations are then ref l ected i n 
total line flux in 0.1 - 12.8 keV energy band (j Jovanovic and Popovid l2008bl ). 

Thus, this perturbing model could explain the variable "red" feature of 
the Fe Ka line observed in NGC 3516. Besides, the realistic durations of disk 
emissivity perturbations are also obtained i f a central supermassive black h ole 
with mass Mbh = 1 x lO^M© is assumed (I Jovanovic and Popovid . l2008bl ). 



^We should note here that we did not take into account the rotation which can produce 
periodical peaks in the light curves, as it can be seen from Fig. [12] 
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6.2 Absorption by warm X-ray absorbers 



The X-ray emission of AGN could be also significantly absorbed by an out- 
flowing wind, especially in case of so-called Low Ionization Broad Absorption 
Line (LoBAL ) quas ars. Recent observ ations of such quasar s (e.g. Mrk 231 
(iBraito et all . l2004h and H 1413+117 dChartas et al.l. 12007111 c onfirmed the 
presence of X-ray absorbers in these objects. Wang et al. ( 200ll ) detected an 
absorption line at 5.8 keV in nearby {z = 0.0033) Seyfert 1.5 galaxy NGC 
4151. A variable ab sorption line at the same energy has been discovered by 
Nandra et al.l (1l999l ) in NGC 3516 (Fig. [T5l) and was interpreted as a Fe K 
resonant absorption line, redshifted either by inf ailing absorbing material or 
by strong gravity in the vicinity of the black hole. 

Done et al.l (120071 ) found an evidence for a P Cygni profile of the Fe 
Ka line (Fig. [TBI) in narrow line Seyfert 1 galaxies. According to these 
authors, complex X-ray spectra of these objects show strong "soft excess" 




Energy (keV) 



Figure 15: Absorption component at 5.8 keV in the red part of the Fe Ka line 
of Se yfert 1 galaxy NGC 3516, observed by ASCA satellite (jNandra et al 



19991 ). The dashed line shows best fit with model of an accretion disk around 
a rotating (Kerr) black hole. 
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below 2 keV and a sharp drop at ~ 7 keV which can be explained either by 
reflection or by absorption from relativistic, partially ionized material close 
to the black hole. They showed that a sharp feature at ~ 7 keV results from 
absorption/scattering/emission of the iron Ka line in the wind. In the case 
of IH 0707-495 (Fig. \T6\i. this abso r ption feature can be satisfactorily fitted 



by the P Cygni profile (jPone et al.l . 120071 ) 




0.5 1 2 5 

channel energy (keV) 



10 



Figure 16: The full XMM-Newton spectrum of IH 0707-495 wit h the best-fit 



mode l which involves a P Cygni profile for the iron features (IDone et al. 
20071 ). The lower panel shows residuals to the fit. 



There are different models of the X-ray absorbing/ obscurin g regions, like 
absorbing medium comprised of cold absorbing clo udlets bv iFuerst &: Wu 



(2004), but here we will focus on the model given by lJovanovic and Popovic 
(120071 ). developed in order to study how much warm absorbers can change 



the Fe Ka spectral line profile, emitted from a relativistic accretion disk. In 
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this model, absorption region is considered to be composed of a number of 
individual spherical absorbing clouds with the same small radii (see Fig. [T7] 
left), scattered in space so that projections of their centers to the observer's 
sky plane (Xj, 1^) have bivariate normal distribution Af2{fi,T,). Here, /i = 



X 



/iy] and 



paxfTY 



pax cry 



a 



Y 



where nx and py are the means of 



Xi and Yi, ax and ay are their standard deviations and p is the correlation 
between them. Absorbing region presented in Fig. [T7| (left) is obtained for 
the following parameters: p = 0, px = ^A, Py = Ya and ax = ay = Ra, 
where {X^, Ya) is the center and Ra is the radius of projection of entire 
absorbing region. 




X (in Rg) 



Energy (keV) 



Figure 17: Left: Relativistic accretion disk in the Schwarzschild metric par- 
tially covered by a cloud of absorbing material (randomly scattered gray 
dots) . Right: Comparison between the unabsorbed Fe Ka spectral line profile 
(solid line) and corresponding absorbed profile (dotte d line) caused by the ab- 
sorbi ng/obscuring region presented in the left panel (jJovanovic and Popovid . 
2007f ). 



The absorption coeffici ent A(X, Y) for every spheric al cloud in the ab- 
sorbing region is given by (jJovanovic and Popovid . 120071 ) : 



AiX,Y) = il-lAiX,Y))-e 



g{X,Y) Eo-Ea 



(24) 
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where absorption intensity coefficient Ia{X, Y) describes the distribution of 
absorption over the whole region, Ea is the central energy of absorption and 
aE is the width of absorption band (velocity dispersion). 

A comparison between the unabsorbed Fe Ka spectral line profile and the 
corresponding absorbed profile obtained by this absorption model is given in 
Fig. [T7| (right). As it can be seen from Fig. [TTl when the X-ray radiation 
from approaching side of the disk is significantly absorbed/obscured by the 
absorbing region, there is a very strong absorption of the iron line. In such 
case the emission Fe Ka line looks redshifted at ~ 5 keV and is followed by a 
strong absorption line at ~ 7 keV (Fig. [T7]right), which indicates the P Cygni 
profile of the iron line. Thus, this model can satisfactorily explain the P Cygni 
profile of the Fe Ka line in the case when approaching side of the accretion 
disk is partially blocked from our view by the X-ray absorbing/obscuring 
material, while the rest of the di s k is le ss absorbed/obscured and therefore 



is visible (jjovanovic and Popovid . 120071 ) . 



6.3 Gravitational microlensing 

Some recent observational and theoretical studies suggest that gravitational 
microlensing can also induce variability in the X-ray emission of AGN, espe- 
cially in the case of gravitationally lensed quasars. The phenomenology of 
gravitational lensing effects and introduction to this field has been given in 
several review papers and b ooks. References to some of them can be found 



in e.g. iPopovic et al.l (120021 ). Therefore, in this section we will not discuss all 
aspects of such effects in the universe, but instead, we will briefly present the 
basic concepts of gravitational microlensing theory, as well as some examples 
for microlensing influence on the X-ray emission from AGN. 

Gravitational lensing is an universal natural phenomenon where the grav- 
itational force of lensing object induce either the amplification of some back- 
ground source (microlensing), or the appearance of its multiple images 
(macrolensing), due to light bending in a gravitational field of the deflec- 
tor. Separation angle between the images of background source depends on 
the mass of gravitational lens and therefore, multiple images can be observed 
only in the case of a massive lensing object, such as galaxy. One of the most 
famous multiple image lens systems is quasar QSO 2237+0305, also known 
as Einstein Cross, which is located at redshift z = 1.695 (see Fig. [T8|) . Its 
four images are due to lensing effect of galaxy ZW2237+030, located be- 
tween us and the quasar at redshift z = 0.0394. In the case of a small mass 
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lens (e.g. a star), the separation angle is also small and therefore, different 
images of background source cannot be resolved. Instead, its intensity is am- 
plified, causing the changes in the observed light curve. Common name for 
both, macrolensing (or simply lensing) and microlensing is strong lensing, 
contrary to weak lensing, which causes distortions in observed images of 
distant objects which can be then used for studying the mass distribution 
along the line of sight, but such phenomena will not be discussed here. 



D 



Figure 18: Gravitationally lensed system QSO 2237+0305 (Ein 
stein Cross) observed by the Advanced CCD Imaging Spectrom 



eter (ACIS) onboard the Chandra X-ray Observatory ( iDai et al. 
2003h . The green circles in each image are the corresponding Huh 



hie Space Telescope (HST) image positions provided by CASTLES 



( http: / / cfawww.harvard.edu/glensdata/Individual / Q2237.html ) . 



Zakharov et al.l (120041 ) found that cosmologically distributed microlenses 
could significantly contribute to the X-ray variability of high-redshifted {z > 
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2) QSOs. Indeed, microlensing of the Fe Ka line has been repo r ted a t 



least in three ma crolensed QSOs: MG J0414+0534 fIChartas etaP 1200^1 



QSO 2237+0305 flDai et all 120031). and H 1413+117 flOshima et al. 



2001 



Chart as et al.l . |200J). IChartas et al.l (120021 ) detected an increase of the Fe 



Ka equivalent width in the image B of MG J0414+0534 which was not fol- 
lowed by the continuum and explained this behavior by assumption that the 
thermal emission region of the disk and the Compton up-scattered emission 
region of the hard X-ray source lie within smaller radii than the iron-line 
reprocessing reg ion. Analyzing the X-ray variability of QSO 2237+0305A, 
Dai et al.l (120031 ) also measured amplification of the Fe Ka line in component 
A of QSO 2237+0305 but not in the continuum and suggested that the larger 
size of the continuum emission region in comparison to the F e Ka emission 
region could explain this result. In the case of H 1413+117, 1 Chart as et al. 
(120041 ) found that the continuum and the Fe Ka line were enhanced by a dif- 
ferent factor. For more detailed discussion about observational and theoreti- 
cal i nvestigations of t he X- r ay variabili t y in g r avitationally lensed Q S Os se e 
e.g. Jovanovic et al. ( 2003 ): Jovanovic ( 2005[l: Jovanovic fc Popoyic (2005) 



Jovanovid(l2006h : iJovanovic et al.l (120081 ) and lPopovic et al.l (l2006al ). and also 
references therein. 

The influence of gravitational microlensing on the X-ray emission from a 
compact accretion disk of AGN can be studied using numerical simulations 
based on ray-tracing method described in §5.3 of this chapter. In the fol- 
lowing text we will briefly review three models of gravitational microlenses 
(point-like, straight-fold caustic and quadruple microlens), as well as their 
influence on the X-ray radiation of AGN. 

If X and Y are the impact parameters (i.e. the coordinates) which 
describe the apparent position of each point of the accretion disk image 
on the celestial sphere as seen by an observer at infinity, then the am- 
plified brightness (due to gravitational micr olensing influence) is given by 
( Jovanovic et al. . 20031 : Popovic et al.l . 2003al lbl): 

I,{X, Y) = e{X, Y)g\X, Y)S{x - g{X, Y))A{X, Y), (25) 

where x = I'obs/^o {^o and Uobs are the transition and observed frequen- 
cies, respectively), g = I'obs/^em is the shift due to relativistic effects {uem is 
the emitted frequency), e is the disk emissivity and A is the amplification 
caused by microlensing. Usually, we do not know the exact form of A, so we 
are forced to consider the different approximations in order to estimate this 
quantity. 
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The first of tliese approximations is called point-like microlens, and 
it is applied when an isolated compact object (e.g. a star) plays a role of 
gravitational microlens. S uch microlens is chara cterized by its Einstein Ring 
Radius in the lens plane ( Jovanovic et al. . 20081 ): 



ERR 



c2 

or by the corresponding projection to the source plane: 



(26) 



^0 



Di 



'-ERR 



AGm DsDis 



D, 



(27) 



where G is the gravitational constant, c is the speed of light, m is the mi- 
crolens mass and Di, Dg and Dig are the cosmological angular distances 
between observer-lens, observer-source and lens-source , respectively. In this 
case, the amplification is given by the following rel ation (INarayan fc Bartelmann 
19991 : IPopovic et aP . I2OOII : Ijovanovic et~aD . l2003h : 



where 



A{X, Y) 



u{X,Y) 



u'{X,Y) + 2 



u{X,Y)^u\X,Y)+A' 



v/(x - Xo)2 + (y 



Fn)2 



^0 



(2^ 



(29) 



corresponds to the angular separation between the microlens and a source. 
Xq, Iq are the coordinates of the microlens with respect to the disk center 
(given in Rg). 



Th e total observed flux is then given by (1 Jovanovic et al.l . 120031 : IPopovic et al. 
2003al Jbl): 



F{x) 



Ip{x)dfl 



(30) 



image 



where dfl is the solid angle subtended by the disk in the observer's sky and 
the integral extends over the whole disk image. 

An illustration of a point-like gravitational microlens crossing over an 
accretion disk in t he Kerr metric with angular momentum a = 0.998 is pre- 
sented in Fig. M flPopovic et all . l2002l . l2003bl l The simulated unamplifled 
and amplifled Fe Ka line proflles due to dif ferent positions of the microlens 
are shown in Fig. [201 (IPopovic et al.l . l2003al ) . 
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Figure 19: Illustration of a point-like gravitational microlens crossing over 
an accretion disk in Kerr metric with angular momentum a = 0.998 
(jPopovic et al.l . 120021 . l2003bl ). Einstein Ring of the microlens is schemati- 



cally presented by yellow Euclidian circle. 



From Fig. [201 it can be seen that point-like microlens could induce strong 
changes of the Fe Ka line shape and intensity, depending on the location of 
the microlens. In the first place, we have in mind the changes in number of 
peaks, their relative separation and the peak velocity. Secondly, such transit 
of the point-like microlens could cause an asymmetrical enhancement of the 
line profile, and the maximum of the amplification would be obtained for 
negative values of x-coordinate, i.e. for the approaching side of the rotating 
accretion disk, and therefore it would affect mainly the blue part of the Fe 
Ka line. 

In most cases we cannot simply consider that microlensing is caused by 
an isolated compact object, but we should take into account that such micro- 
deflector is located in an extended object (typically, the lens galaxy). There- 
fore, when the size of the microlens projected Einstein Ring Radius is larger 
than the size of the accretion disk and when a number of microlenses form 
a caustic net, we de scribe the micro l ensing effects by a straight -fold caus- 
tic approximation (jjovanovic et al.l . l2008l : iPopovic et al.l . l2003al ). where the 
amplification at a point of an extended sou rce (accretion disk) close to the 
caustic is given by (IChang fc Refsdall . Il984l ): 



A{X, Y)=Ao + K 



^caustic 



(31) 



In above expression Aq is the amplification outside the caustic and K = A^jS 
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Figure 20: The unamplified (dotted line) and amplified Fe Ka line profiles 
(solid line) for different positions of the point-like microlens, which coordi- 
nates (expressed in R n) in respect to the a ccretion disk center are given at 



the top of each panel (jPopovic et al.l . l2003al ). The simulations are performed 
for an accretion disk in Kerr metric with angular momentum a = 0.998 and 
with the following disk parameters: Rj„ = 6 Rg, Rout = 20 Rg, i = 35°, 
q=2.5. The Einstein Ring Radius of the point-like microlens is 10 R^. The 
relative intensity (y-axis) is in the range from to 3, and energy shift g due 
to relativistic effects (x-axis) is in the range from 0.4 to 1.2. 
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is the caustic amplification factor, wliere (3 is constant of order of unity. ^ 
is the distance perpendicular to the caustic and C,c is the minimum distance 
from the disk center to the caustic (both in Rg). The "caustic size" r caustic is 
the distance ^ for which the caustic amplification is equal to 1, and therefore 
this parameter defines a typical linear scale for the caustic. H[k{S, — ic)) is 
the Heaviside function which equals 1 for — ^c) > and otherwise it is 0. 
K is ±1 depending on the direction of caustic motion: if the direction of the 
caustic motion is from approaching side of the disk k = — 1, otherwise it is 
+1. Also, in the special case of caustic crossing perpendicular to the rotating 
axis n = +1 for direction of caustic motion from —Y to +Y, otherwise it is 
-1. 

A microlensing event where a c austic crosses oyer an accretion disk can 
be described in the following way fjjovanovic et al.l . l2008l ): before the caus- 
tic reaches the disk, the amplification is equal to Aq because the Heaviside 
function of Eq. ( 13T1) is zero. Just as the caustic begins to cross the disk, 
the amplification rises rapidly and then decays gradually towards Aq as the 
source moves away from the caustic-fold. 

An illustration of the straight-fold caustic crossing over an accretion disk 
in Schwarzschild metric is presented in Fig. [211 and the corresponding effects 
on t he shapes of the X-ray continuu r n and t he Fe Ka line are shown in Fig. 
(jJovanovig . l2006l : iPopovic et al.l . l2003bl ). As it can be seen from Fig. 
caustics could also induce significant amplifications of both, the X-ray 
continuum and Fe Ka line. 




Figure 21: Illustration of a straight-fold caustic gravit ational mi c rolens 
crossing over an accr etion disk in Schwarzschild metric (jJovanovig . 12006 : 
Popovic et aD . l2003bh . 
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Figure 22: The deformations of the X-ray continuum (left) and the Fe Ka 
line (right) due to influence of straight-fold caustic from Fig. [21] Undeforme d 
profiles are presented by solid and deformed by dotted lines ( Jovanovi^ . I2OO6I ) . 



The most complex approximation for gravitational microlensing amplifi- 
cation is so called quadruple microlens. This model is applied to obtain 
a spatial distribution of magnifications in the source plane (where an accre- 
tion disk of AGN is located), produced by a random star field placed in the 
lens plane. Such spatial distribution of magnifications is called microlens- 
ing map, microlensing pattern or caustic network. If we consider a 
set of N compact objects (e.g. stars) which are characterized by their posi- 
ti ons Xj and their m asses rrii, then the normalized lens equation is given by 



flAbajas et all . 120071 ): 



N 

E 

i=l 



rrii 



+ 



1 - + 7 

1 - - 7 



(32) 



where x and y are the normalized image and source positions, respectively. 
The sum describes the light deflection by the stars and the last term is a 
quadruple contribution from the galaxy containing the stars, where 
smooth surface mass density and 7 is an external shear. The total surface 
mass density or convergence can be written as k = + k^, where k^, rep- 
resents the contribution from the compact microlenses. The corresponding 
microlensing map is then deflned by two parameters: the convergence - k, 
and the shear due to the external mass - 7. 

For some speciflc lensing event one can model the corresponding magnifl- 
cation map using numerical simulations based on ray-shooting techniques (see 
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e.g. Jovanovic et al. ( 2008 ): Popovic et al. f 2006al Jbl) and references therein), 
in which the rays are shot from the observer to the source, through the ran- 
domly generated star field in the lens plane. An example of a magnification 
map for a "typical" lens system (i.e. for a lens system where the redshifts of 
microlens and source are zi = 0.5 and Zg = 2, respectively) is presented in the 
left panel of Fig. [231 and the corresponding X-ray, UV and o ptical continuum 



variations are given in the right panel of the same figure (iJovanovic et al. 



20081 ). The size of this magnification map is IQ rjo x IQ rjo and it is generated 
using the following arbitrary values for convergence and shear: n = 0.45 and 
7 = 0.3. The light curves in the right panel of Fig. [23] are produced when an 
accretion disk crosses over the magnification pattern in the left panel, along 
the path denoted by the white solid line. 





-looao 10000 

Distance from Ihe center of accretion disc {in R J 



Figure 23: Left: Magnification map of a "typical" lens system, where the 
white solid line represents a path of an accretion disk center. Right: Varia- 
tions in the X-ray (solid), UV (dashed) and optical (dotted) spectral bands 
corresponding to the path in the magnification map of a " typical" lens system 
(left). 



As one can see from Fig. [23] (right), the variations of the X-ray continuum 
due to microlensing are much stronger and faster in comparison to the varia- 



tions i n UV and optical spectral bands (IJovanovic et al 
2006al Jbll Usually, so called High Amplifi cation Events, i.e 



20081 : ]Popovic et al. 



peaks in the light curves are also analyzed (iJovanovic et al. 



asymmetric 
20081 ). The re- 
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suits presented in Fig. [23] (right) show that the rise times of such events are 
the shortest and their frequency the highest in the X-rays, in comparison to 
the UV/optical spectral bands. 

In or der to explain the obse rved X-ray variability in gravitationally lensed 
quasars, iPopovic et al.l ( j2006al ) considered a microlensing magnification map 
for a specific case of the QSO 2237+0305A image (Fig. [2ll top). This map 
with 1 ?7o X 2 ?7o on a side was generated using the following parameters: 
K = 0.36, 7 = 0.40 and the mean mass of randomly distributed deflectors 
(m) = 0.35 Mq. The corresponding variations of the total Fe Ka line and 
X-ray continuum fluxes are presented in Fig. [21] (bottom). 

As one can see from Fig. [2l] (bottom), there is a global correlation be- 
tween the total line and continuum fluxes during the complete simulated 
microlensing event, but while the continuum flux v ariations a r e smo oth and 
slow, the line flux varies much stronger and faster (jjovanovig . [2005]). More- 
over, there is a certain period in the middle of the microlensing event in 
which the Fe Ka line flux changes rapidly, while the continuum flux remains 
nearly constant. This indicates that different behavior of the line and X-ray 



conti nuum variability, observed in some lensed quasars (1 Chart as et al.l . [2002 



20041 ). may be explained by microlensing hypothesis if the line originates 



from the innermost part of the accretion disk, while the X - ray continuum is 



emitted from a larger region in the disk ([Jovaiiovic et al.[ . [2003[ : [Jovanovid . 



20051 : [Jovanovic fc Popovid . [20051 : [Popovic et al.[ . [2006a[ ) . 



On the basis of all facts presented in this section, one can conclude that 
gravitational microlensing can produce signiflcant variations and ampliflca- 
tions of the Fe Ka line and X-ray continuum. During a microlensing event, 
even very small mass objects could produce noticeable changes in the X- 
ray radiation from accretion disks of AGN. These changes are signiflcantly 
larger than the corresponding effects in the optical and UV emission. Also, 
microlensing hypothesis can satisfactorily explain the excess in the iron line 
emission, observed in some gravitational lens systems. 



7 Conclusion 

In this chapter we discussed the X-ray radiation from relativistic accretion 
disks around supermassive black holes, supposed to exist in the centers of all 
Active Galactic Nuclei. The X-ray radiation is created when the accreting 
material loses its angular momentum due to friction between the adjacent 
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Figure 24: Top: Microlensing magnification map of QSO 2237+0305A image 
witli 1 ?7o X 2 ?7o on a side and a sclieme of an accretion disk with separated Fe 
Ka line and X-ray c ontinuum emit t ing reg ions which outer radii are 20 and 



100 Rg, respectively (jPopovic et al.l . l2006al ). The straight horizontal line rep 



resents a path of the disk center. Bottom: The corresponding amplifications 
of the Fe Ka line and X-ray continuum total fluxes for different positions of 
the disk center along the path on the magnification map in the left panel. 



layers in the disk and spiral towards the central black hole, releasing the grav- 
itational energy. Part of this energy increases the kinetic energy of rotation 
and the rest is turned into the radiation, which is then emitted from the disk 
surface in the X-ray, UV and optical spectral bands. It is the only known 
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mechanism of converting gravitational potential energy into radiation, which 
is sufficiently efficient to explain the high luminosities of some observed AGN. 

The standard model of an accretion disk, including its emission, accretion 
rate, luminosity, structure and spectral distribution was described in details. 
A significant part of this chapter was devoted to modeling of the observed 
X-ray radiation in both the Fe Ka spectral line and X-ray continuum. The 
modeled X-ray radiation is then used to analyze the several effects of strong 
gravity in the vicinity of black holes which have the influence on the shapes 
and intensities of the observed Fe Ka line and X-ray continum. We showed 
that the X-ray emission strongly depends on the angular momentum and 
mass of the central black hole, as well as on several parameters describing 
the X-ray emitting region in the disk. This was demonstrated on several 
examples, showing various aspects of the X-ray emission from the accretion 
disk around Schwarzschild and Kerr black holes. Since the X-ray radiation 
is probably produced in a very compact region near the black holes, such 
comparison between observed and modeled X-ray radiation can provide us 
some essential information about the plasma conditions and the space-time 
geometry in the vicinity of those black holes. 

We also studied some possible causes for the rapid and irregular variability 
of the Fe Ka line and X-ray continuum, observed in the X-ray spectra of 
some AGN. Our results show that, at least, three mechanisms could explain 
such variability: perturbations of disk emissivity, absorption by warm X-ray 
absorbers and gravitational microlensing. 

Taking all this into account, we can conclude that the X-ray emission 
from accretion disks of AGN can be used as a powerful tool for reveahng 
the physics and geometry in the vicinity of their central supermassive black 
holes. But, to use the Fe Ka line profile for such investigations, one should 
also take into account the effects which can disturb its profile (such as e.g. 
absorption, perturbations in the disk and gravitational lensing). 
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